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Abstract: Realizing a quantum theory for gravity based on Asymptotic Safety hinges 
on the existence of a non-Gaussian fixed point of the theory’s renormalization group flow. 
In this work, we use the functional renormalization group equation for the effective average 
action to study the fixed point underlying Quantum Einstein Gravity at the functional 
level including an infinite number of scale-dependent coupling constants. We formulate a 
list of guiding principles underlying the construction of a partial differential equation en¬ 
coding the scale-dependence of /(i2)-gravity. We show that this equation admits a unique, 
globally well-defined fixed functional describing the non-Gaussian fixed point at the level 
of functions of the scalar curvature. This solution is constructed explicitly via a numerical 
double-shooting method. In the UV, this solution is in good agreement with results from 
polynomial expansions including a finite number of coupling constants, while it scales pro¬ 
portional to R^, dressed up with non-analytic terms, in the IR. We demonstrate that its 
structure is mainly governed by the conformal sector of the flow equation. The relation of 
our work to previous, partial constructions of similar scaling solutions is discussed. 
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1 Introdnction 

Renormalization group (RG) flxed points are crucial for determining the properties of the 
underlying quantum held theory. They are lighthouses on theory space because in their 
vicinity one can explore the corresponding quantum held theory via scaling techniques 
and, eventually, via conformal held theory methods. Fixed points also play a crucial role 
in providing a general notion of renormalizability and ultraviolet completion of a quantum 
theory. In fact, in the vicinity of an ultraviolet (UV) flxed point the RG flow can be 
linearized and characterized by directions along which it is attracted to or repelled from 
the flxed point for increasing energies, thus providing a generalized notion of relevant and 
irrelevant deformations of the theory. The UV completeness of the theory is then ensured 
by finding a given RG trajectory that is attracted towards the flxed point at high energies, 
implying that any UV cutoff can be removed in such a way that all essential dimensionless 
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coupling constants remain finite. Such an UV complete theory is also predictive if the 
dimension of the critical hypersurface that is attracted towards the fixed point in the 
UV is finite. In this case, only a finite number of experiments is required to pinpoint 
the specific RG trajectory that completes gravity in the deep UV regime. Depending on 
whether the fixed point corresponds to a free (Gaussian fixed point) or interacting theory 
(non-Gaussian fixed point), the RG trajectories terminating at the fixed point in the UV 
are termed asymptotically free or asymptotically safe, respectively. While it is well known 
that the Gaussian hxed point (GFP) of gravity does not provide a suitable mechanism 
to renormalize the theory, the possibility that the theory possesses at least one suitable 
non-Gaussian fixed point (NGFP) is not speculative and has been explored under various 
convincing approximations. In this work we investigate the Asymptotic Safety mechanism 
for gravity relaxing some of these approximations. In particular we construct the first 
consistent and globally well-defined fixed functional for Quantum Einstein Gravity in four 
spacetime dimensions. 


1.1 Quantum Gravity via Asymptotic Safety 

The RG flow of any theory admits a GFP at which the theory itself becomes non-interacting. 
In the neighborhood of the Gaussian point it is thus possible to use the methods of pertur¬ 
bation theory to investigate the UV completion of the theory under consideration. These 
methods have been applied very successfully to Yang-Mills theories thanks to the property 
of asymptotic freedom. However, it has been known for many years that the GFP cannot 
provide a UV completion mechanism of General Relativity [1-3]. The question whether 
gravity can be formulated as a consistent and predictive (asymptotically safe) quantum 
field theory then becomes closely intertwined with the question whether there are other 
fixed points of the RG on the gravitational theory space which could provide such a UV 
completion. Indeed it has been conjectured already in the late seventies that there exists a 
suitable NGFP on the gravitational theory space which could render gravity asymptotically 
safe [4]. 

The primary tool for investigating the Asymptotic Safety mechanism has been a 
Wilsonian-type functional renormalization group equation (FRGE) [5-7] 


kdkTk = ^Tr 




-1 

kdkl^k j 


( 1 . 1 ) 


( 2 ) 

in which F^ denotes an effective “average” action, F^ ' is its second variation with respect 
to the fields, and T^k is a suitable infrared-regulator which ensures that the rhs of the 
equation is finite and peaked around the renormalization group scale k. The effective 
average action T k represents an effective action of the system in which the UV modes have 
been integrated out in the path-integral using the scale k as a reference, thus providing 
a natural effective action for processes occurring at energies E k. In the limit fc —>• 0 
the IR cutoff IZk vanishes and the effective average action coincides with the full effective 
action of the system F = Ffc=o- The flow (1.1) is exact in the sense that no approximation 
is used in its derivation. 
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Notably, the FRGE allows to compute approximate RG flows of a theory without the 
need of a small expansion parameter. For gravity this approach has been pioneered in [8] 
which, by now, has evolved in a very successful research program, see [9-12] for reviews. 
The key strategy for obtaining non-perturbative information on the RG flow is to project 
the exact flow on a subspace spanned by the “most relevant interactions” for the physics at 
hand. Practically, this strategy is implemented by making a suitable ansatz for F^ limiting 
the set of interaction terms to those which are argued to describe the relevant physical 
information. The beta functions, which describe the scale-dependence of the coupling 
constants contained in the ansatz and drive the approximated RG flow, are constructed 
by substituting the ansatz for F^ into the FRGE and retaining the terms contained in the 
ansatz itself only. While these beta functions lack the control parameter h of the loop 
expansion of perturbation theory [13], they are fully non-perturbative in the couplings and 
thus, for sufficiently big truncations, reliable also away from the GFP. 

The projection of F^ to finite dimensional subspaces has led to a substantial body of 
evidence supporting the existence of a gravitational NGFP suitable for Asymptotic Safety. 
Starting from the Einstein-Hilbert truncation, which contains scale-dependent Newton’s 
and cosmological constants only [8, 14, 15], the existence of the NGFP has successively 
been established on truncations of increasing complexity and field content, see [16-19] for 
selected references. For the case in which the gravitational degrees of freedom are carried by 
a (Euclidean) spacetime metric, this line of research has already established striking results. 
Foremost, all computations corroborated the existence of a NGFP with suitable properties 
for Asymptotic Safety. In particular, the NGFP has been put to test and proved robust 
to the inclusion of the problematic operators appearing in the perturbative counter-terms 
in Ffc [20, 21], thus directly showing the feasibility of the program beyond perturbation 
theory. Despite this fixed point being associated with an interacting quantum field theory in 
which the (background) Newton’s constant scales with an anomalous dimension = —2, 
classical power-counting may still constitute a good ordering principle for identifying the 
relevant deformations of the hxed point as shown in [22-25]. In particular the number of 
relevant parameters of the theory could be as low as three, implying that, compared to 
many other quantum field theory models, an asymptotically safe theory of gravity can be 
very predictive. Moreover, stringent bounds on the numbers of free matter fields compatible 
with the Asymptotic Safety mechanism have been obtained in Refs. [26, 27]. At the level 
of the formalism, there has also been progress on computing properties and expectation 
values of the fluctuation fields around a given background [28-32] . By now it is appreciated 
that this type of computations, known as bi-metric computations, may be essential for 
precision computations of the properties of NGFP such as its critical exponents, addressing 
global aspects of gravitational RG flow by, e.g., establishing a suitable c-theorem [33], or 
establishing background-independence of the formalism [34] . 

While the case where the gravitational degrees of freedom are carried by metric vari¬ 
ables is understood the best, this scenario is not the only way to realize Asymptotic Safety 
in gravity: a suitable NGFP has been found to occur rather independently of the type 
of gravitational degrees of freedom. Investigations of the gravitational RG flows based 
on the ADM-formalism [35, 36], the Einstein-Gartan formalism [37-40] and the “tetrad 
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only” formalism [41, 42] also identify a NGFP whose properties are strikingly similar to 
the one encountered in the metric formalism. This is rather remarkable since, while these 
formulations are equivalent at the classical level, it is far from clear that this equivalence 
also holds at the quantum level. Moreover, there are also first indications that Asymptotic 
Safety may be realized in unimodular gravity [43-45]. This leads to the perhaps surprising 
conclusion that, at the level of finite-dimensional projections, the appearance of a NGFP 
in gravitational theories is rather commonplace instead of exceptional.^ 

1.2 From fixed points to fixed functions 

Currently, one of the key conceptual challenges in the Asymptotic Safety program is the 
leap from approximations for F^ containing a finite number of scale-dependent couplings 
to truncations containing an infinite number of coupling constants. In the latter case, F^ 
contains scale-dependent functions. Substituting such an ansatz into the exact functional 
renormalization group equation leads to a (system of) partial differential equations (PDF) 
which govern the scale-dependence of these functions. The fixed points appearing in finite- 
dimesional approximations of F^ are then promoted to fixed functions, namely global k- 
stationary solutions of the PDF. (Since both fixed points and fixed functions induce fixed 
functionals, we use the notion of fixed functions in order to stress that we work in the 
realm of infinitely many scale-dependent coupling constants.) 

The simplest setting where these generalizations can be implemented for gravity in¬ 
cludes an arbitrary function of the scalar curvature R 

rr''[9]= / d^x^gUR), (1.2) 

where g denotes the (Fuclidean) spacetime metric. This so-called /(i?)-truncation provides 
the simplest example for a truncation including a scale-dependent function, since the PDF 
governing the ^-dependence of fkiR) can be derived by using a maximally symmetric 
background with covariantly constant curvature. 

Starting from the flow equation derived in [23], the existence of fixed functions f*{R) 
have been investigated by various groups in d = 3 [52-54] and d = A [55-59] spacetime 
dimensions. Quite unsettling, the verification of a suitable NGFP at the level of fixed func¬ 
tions turned out to be extremely challenging: while the finite-dimensional computations 
always produced a suitable UV fixed point regardless of the computational details and 
setting, up to now the fixed function completing the four-dimensional NGFP in the ansatz 
(1.2) is still elusive. This deficit acted as catalyst to readdress some of the fundamental 
questions in the functional renormalization group approach to Quantum Finstein Gravity 
(QFG). 

Substituting the ansatz (1.2) into the FRGF and projecting on the corresponding 
subspace leads to a non-linear PDF which is first order in k and third order in the R- 
derivatives. Imposing the condition that the fixed function is /^-stationary, this equation 

^Notably, Hofava-Lifshitz gravity [46] also possesses a UV fixed point, which is different from the one 
underlying the Asymptotic Safety scenario [47]. It is conceivable that the Monte Carlo simulations carried 
out within the Causal Dynamical Triangulations Program [48-51] are capable of probing both the Hofava- 
Lifshitz and Asymptotic Safety universality class. 
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reduces to a non-linear ordinary differential equation (ODE) of third order. Thus, locally, 
there is a three-parameter family of solutions to this equation. It was then noticed early- 
on that the condition for the solution to exist globally may reduce the number of free 
parameters and may produce isolated fixed functions if there is a balance between the 
order of the ODE and its singularity structure [56] 

degree of ODE — number of singularities = 0 . (1-3) 

Thus any singularity reduces the number of free parameters by one. While the number of 
fixed singularities can essentially be read off from the ODE, determining the occurrence of 
movable singularities often requires studying explicit solutions. Thus, so far, only fixed sin¬ 
gularities played a role in restricting the number of free parameters in the fixed function in 
the /(i?)-truncation. As a important consequence of (1.3) it was then understood that the 
initial flow equation contains too many fixed singularities in order to admit global solutions 
[55, 56]. Based both on physical arguments and on the experience maturated during the 
course of several /(i?)-type analysis, we propose the following guiding principles towards 
the construction of a suitable hxed function in four-dimensional QEG, which will turn out 
to be a set of sufficient conditions for the existence of said fixed function: 

1) The background should admit only one topology. 

The initial flow equation has been derived on a maximally symmetric, compact background 
sphere S'^ with positive scalar curvature. A detailed comparison of the flow equations for 
fk{R) in the three-dimensional conformally reduced setting established that the topology of 
background (non-compact with negative scalar curvature or compact with positive scalar 
curvature) manifestly influences the singularity structure of the PDE. Additionally, the 
flow equation on a hyperbolic background with negative scalar curvature does not follow 
from the analytic continuation of the flow equation on a sphere. We thus argue that it is 
potentially inconsistent to include more than one topological sector in the flow equation 
avoiding topology changes in the background. In the sequel, we restrict the construction 
to the domain R > 0 and do not demand that the solutions admit an continuation to the 
R < 0 region, thereby avoiding that the fix functional has a suitable continuation into a 
domain with different topology. 

2) The path-integral measure should not introduce spurious singularities. 

As it turns out, the measure in the gravitational path integral has a crucial effect on the 
singularity structure of the ODE. Using a linear background held decomposition 

Qfj.i' ~ 9iJ,u (f'4) 

together with a standard measure for the metric huctuations and evaluating the how 
equation based on the transverse-traceless decomposition of generally introduces hxed 
singularities from the contributions of the gauge degrees of freedom, ghost terms, and Ja- 
cobians arising in the held decomposition. We believe that these singularities have to be 
considered as spurious, since their number can be changed at will by varying unphysical 
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parameters. We will resort to the flow equation derived in [60] in order to avoid the spuri¬ 
ous poles of the gauge-fixing sector and of the Jacobians of the held decompositions. This 
framework can be motivated by geometrical considerations, since it corresponds to a trivi- 
alization of the gauge-bundle in held space, and is equivalent to a particular simultaneous 
choice of regulator function TZk and of gauge-hxing, which allows for a precise cancellation 
between the functional traces capturing the gauge-degrees of freedom and ghost contribu¬ 
tions. 

3) The asymptotic behavior is captured only by the exact heat-kernel. 

The rhs of the FRGE (1.1) is a trace involving a (generally complicated) differential op¬ 
erator. In hnite-dimensional truncations, this functional trace is typically evaluated by 
applying the early-time expansion of the corresponding heat-kernel. This is the natural 
approach if T^ is organized via a derivative expansion, for example in powers of the back¬ 
ground curvatures. It turns out, however, that at the level of full functional truncations 
the use of the early-time expansion and the evaluation of the trace based on the exact 
heat-kernel on leads to a different asymptotic behavior when <C R. This difference 
is caused by non-analytic contributions to the heat-kernel, related to the effect that on 
the diffusing particle may return to its origin multiple times. This effect is not captured by 
the (analytic) early-time expansion. These terms are, however, essential for reproducing 
the correct late-time behavior of the heat-kernel (and are yet another manifestation of the 
compactness of the background manifold). Therefore, we base the evaluation of the oper¬ 
ator trace (1.1) on the exact heat-kernel including non-analytic contributions to correctly 
account its asymptotic behavior. 

4) Scalar and TT-fluctuations should be integrated out simultaneously. 

The initial flow equation has been derived using the Laplacian of the background sphere 
for the construction of a covariant regulator. In the nomenclature of [24], this corresponds 
to a type I cutoff. Our flow equation generalizes this construction by allowing a non-zero 
endomorphism in the regulator function, using a so-called type II cutoff. Besides making 
the solutions of the ODE more stable, the freedom of the endomorphism allows to satisfy 
both the singularity count (1.3) and the correct asymptotic behavior of the solution. More¬ 
over, choosing different endomorphisms in the different sectors of the flow equation allows 
adjusting the scale at which the lowest eigenmode of the fluctuation field is integrated 
out. We thus need a physically motivated principle to fix this additional freedom. We 
will impose the principle of equal lowest eigenvalues [60] , enforcing that the lowest energy 
mode in the transverse-traceless and scalar sectors are integrated out at the same value of 
k. Since the order of the highest derivative of fk{R) differs in the transverse-traceless and 
scalar sector of ODE, this choice in particular ensures that the order of the ODE does not 
change on the integration interval. 

5) All quantum fluctuations need to integrated out. 

A crucial question related to the existence of globally well-dehned solutions is the integra¬ 
tion range for which the solution ought to exist. The ODE describing the hxed function 


- 6 - 


is typically formulated in terms of the dimensionless curvature r = R/k'^. For fixed (back¬ 
ground) curvature R the interval k G [0,oo[ is then mapped to r G [oo,0], so it seems 
natural to demand that solutions are regular on the entire positive real axis. This is in 
agreement with the intuition from scalar field theory in a flat background and the require¬ 
ment that the regulator TZ^ oc vanishes for —)• 0, so that all quantum fluctuations are 
integrated out. Based on this logic the domain of definition for the fixed function may be 
inferred by imposing that 

a) all quantum fluctuations are integrated out, 

b) there are no residual regulator effects in the flow equation. 

Depending on the particular choice of regulator, and in particular for Type II cutoffs in 
combination with the Litim regulator [61], these requirements may already be satisfied 
at a hnite rmax- In particular this construction removes the regulator effects not via the 
“canonical mechanism” sending A; —)• 0 but by the stepfunction containing the eigenvalues 
of the differential operator becoming zero once the last fluctuation mode is integrated out. 
In such a case, after the last eigenfluctuation is integrated out, the flow will be driven by 
the canonical scaling of the function f{R). Based on this mechanism the r-interval, on 
which the global solution is required to be regular, may be restricted to a finite interval 
r G [0, Tmax]) where r max depends on the choice of endomorphism included in the regulator 
and on the regulator itself. If the flow equation is derived using a smooth regulator or 
the behavior of a smooth regulator is mimicked by applying a smoothing procedure the 
contribution of the quantum fluctuations are smeared over the entire half-axis r G [0,oo]. 
Since we apply the smoothing procedure [55] in the construction of our flow equation, we 
consequently require that a proper hxed function is globally well-defined on this domain. 

In this work, we establish that these guiding principles lead to a PDE governing the 
scale-dependence of the ansatz (1.2) which gives rise to an isolated, unique and globally 
well-defined fixed function, thus reconciling previous results based on finite-dimensional 
polynomial truncations with a functional ansatz. While it could have been argued via 
the counting theorem that only a finite number of solutions was present, it is in fact 
amazing to discover that the solution resulting from the application of the above principles 
is unique. Let us stress again, however, that our guiding principles constitute a set of 
sufficient conditions for the existence of said fixed function. At this stage it is unclear 
whether all principles are strictly necessary for establishing a suitable fixed function. In 
particular, the list does not incorporate the recent observation [59, 62] that the combination 
of a flow equation with the equation for the modified split Ward identities may lead to 
substantial simplifications of the PDFs prescribing the fixed functions of the theory. 

The remaining paper is organized as follows. We start by illustrating the effects of 
the singularity theorem (1.3) together with the analytic and numerical techniques used to 
analyze the existence of fixed functions in a simple toy model in Sect. 2. Implementing 
the guiding principles discussed above, the ODE governing the existence of hxed functions 
in the /(i?)-truncation is constructed in Sect. 3. The properties of the resulting ODE 


- 7 - 


are discussed in Sect. 4 and a numerical study of its solutions is carried out in Sect. 5. 
The central result is the isolated, globally well-defined fixed function shown in Fig. 4. We 
end with a discussion in Sect. 6. The analogue of the full fixed function arising from the 
conformally reduced approximation of our flow equation is constructed in Appendix A. 

2 Isolated global solutions: a toy model 

Before embarking on the quest of constructing our flow equation for /(i2)-gravity and its 
fixed functions, it is useful to illustrate the techniques used in the analysis by applying 
them to a simple, completely tractable example. For this purpose, we analyze the second 
order ODE 

/(x) = -^^, y(0) = 7/0 , y'(0) = 2/1 , (2.1) 

on the interval x G [0, oo[. Locally, at a: = 0, the equation admits a two-parameter family 
of solutions which we parameterize by ( 2 / 0 ) 2 /i)- 

Inspecting the rhs of eq. (2.1) one easily sees that the denominator has a root at 
iCsing = 1 which constitutes a fixed singularity. A solution with generic initial conditions 
( 2 / 0 , 2 / 1 ) imposed at x = 0 will have a diverging second derivative at Xging. A necessary 
condition for regularity of the rhs is 


2/(a^)L=x,„, = 0, (2.2) 

i.e., a globally well-defined y{x) has at least a simple root at Xging- Thus it is expected that 
not all solutions can be extended to global solutions which are well-defined on the entire 
positive half-axis. The fixed singularity provides a non-trivial boundary condition which 
puts additional constraints on the parameters (2/0)2/i) characterizing the local solutions. 
We now discuss various techniques which allow to analyze this constraint. 

The exact analytic solution 

The linear ODE (2.1) is still simple enough so that its general solution can be given in 
terms of modified Bessel functions 

y{x) = \/l- X (ci Ii{2y/l - x) + C 2 Ki{2y/l - x)) . (2.3) 

Evaluating the solution in the limit x —>■ 1“, the term containing Ji vanishes while the 
Ki term remains finite. Thus the condition (2.2) implies C 2 = 0. Indeed one can explicitly 
verify that only solutions with C 2 = 0 remain real on x G [0,oo[: if C 2 7^ 0 the iLi-term 
introduces a non-zero imaginary part if x > 1. This property is illustrated in the left 
diagram of Fig. 1. The regularity property (2.2) is thus a sufficient condition for the 
existence of a global real-valued solution in the positive real half-axis. 

Expanding (2.3) in a Taylor-series at x = 0 up to linear order gives the relation between 
the integration constants Cj to the parameters {yo,yi) 

2/o = ci/i(2)+C2iLi(2), 2/i = -ci(/i(2)+/2(2))-C2(iLi(2)-iL2(2)) . (2.4) 
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Figure 1. Illustration of a onc-parameter family of solutions (2.3) given by ci = 1 and various 
values of C 2 (left diagram). The solution C 2 = 0 can be extended beyond the fixed singularity 
s^sing = 1- For C 2 ^ 0 the solutions terminate at the hxed singularity. The right diagram shows 
the values {yo,yi) satisfying the condition (2.2). The bold blue line is obtained from the exact 
solution, the red squares indicate the results from the numerical shooting method, and the green 
line approximates the exact solution by a polynomial with N = 12. 


Imposing the condition for a regular, global solution, C 2 = 0 and eliminating ci from 
the system (2.4) shows that the global solutions are given by a line in the (yo; yi)-pla'ne, 
satisfying 

yi = -2/0 (l + 11 ^) ~ -1.433 yo . (2.5) 

This line is displayed as the bold blue line in the right panel of Fig. 1. Thus the interplay 
between demanding that a solution to exist globally and the fixed singularity of the ODE 
reduces the number of free parameters by one. 


Polynomial expansion at x = 0 

Concerning the ODE describing hxed functions in quantum gravity, it is not likely that 
exact, analytic solutions can be found. Thus other techniques to derive the condition for a 
global solution (2.5) are required. A tractable method approximates the exact solution by 
a polynomial of hnite order 

N 

y{x) ■ ( 2 . 6 ) 

71=0 

Substituting this ansatz into (2.1) and matching the coefficients multiplying powers of x 
up to x^~‘^ yields — 1 equations relating the coefficients 


2/2 = ^ 2/0 
2 / 3=1 (2/0 + 2/1) 


VN 


1 

N{N - 1) 


N-2 

YI y^- 

n=0 


(2.7) 
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This system can be used in two ways. Firstly, it can be complemented by the two 
equations appearing at orders and : 


N-l 

= 0 > 

n=0 


N 


'^yn = 0- 

n=0 


( 2 . 8 ) 


Combining eqs. (2.7) and (2.8) then shows that this completion scheme only leads to the 
trivial solution = 0, V n G 0,..., A^. 

Alternatively, one can follow the strategy of fixing free coefficients by demanding the 
global existence of the solution. For this case, it is convenient to cast the system (2.7) into 
the form 


n — 2 

Un — Un—l T 

n 


1 


.yn- 2 , n>2 (2.9) 

n[n — 1) 

from which the Un can be determined as functions of the parameters yo,yi in a recursive 
manner. Fixing = 12 and evaluating the regularity condition (2.2) for the polynomial 
(2.6) then yields the regularity condition 


1004933203 

701216922 


-1.433yo. 


( 2 . 10 ) 


This condition is illustrated by the thin green line in the right panel of Fig. 1 and is essen¬ 
tially identical to (2.5). Notably, fixing all boundary conditions at x = 0 (first strategy) 
or demanding the existence of a globally well-defined solution (second strategy) does not 
necessarily lead to the same result. 


Numerical shooting method 

The numerical shooting method approximates the space of initial conditions spanned by 
yO) yi by a lattice with a suitable lattice spacing. For illustrative purposes we use Ayi = 0.1, 
corresponding to 4941 pairs of initial conditions on the (yo; ?/i)-plcme shown in Fig. 1. At 
these lattice points, the ODE (2.1) is then solved numerically on the interval x G [0, Xsmg —e] 
employing a standard ODE-solver. Practically, we chose e = 10“^. 

The numerical solutions provide a map 

(yo,yi) ^ e(x;yo,yi)l^=^^_, . (2.11) 

The regularity condition (2.2) corresponds to e(xsing — e, yo,yi) = 0 and defines a submer¬ 
sion on the space of initial conditions (yo,yi)- Technically, this condition is implemented 
by looking for sign changes in e(yo,yi) along the strips of constant yo. The resulting (in¬ 
terpolated) values for yi indicating a zero of (2.11) are given by the red squares displayed 
in the right panel of Fig. 1. Fitting a linear function to the data points yields 

yi«-1.434yo. (2.12) 

This matches the exact result (2.5) with three digit precision. Thus all three methods, exact 
solutions, polynomial approximations and the numerical shooting method are suitable to 
quantitatively analyze the constraints implied by the existence of a global solution on the 
free parameters of the ODE. 
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3 The flow equation for fk{R) 


In this section, we construct the PDE encoding the scale-dependence of fk{R) based on the 
guiding principles laid down in the introduction. The key results are eqs. (3.21) and (3.24) 
which constitute the starting points in our search for fixed functions in the subsequent 
sections. 


3.1 Projecting the RG flow 

The construction of the flow equation (1.1) employs the background field method, splitting 
the physical metric into a fixed background and a fluctuation field 

~ 9iiu T . (3-1) 


When projecting the RG flow onto the ansatz (1.2), it is convenient to choose as 
the metric on the compact four-sphere with background curvature R and preform a 
transverse-traceless decomposition of the metric fluctuations with respect to this back¬ 
ground 

hf_iu — -|- D -|- 2Df^Dy(j -^g^yD cr -|- ^ h. (3-2) 

The component fields consist of a transverse-traceless tensor KJ^y (spin 2), a transverse 
vector , and the two scalars cj, h (spin 0) satisfying the constraints 

D^hly = 0, = = r^h^u = h. (3.3) 


Following [60], the projected flow for the ansatz (1.2) takes the form 


dt^k[9] — 2 ^( 2 ) + 2 ^( 0 ) 

with the explicit expressions for the traces given by 


(3.4) 


T( 2 ) = 6 Tr 
r(o) = 16 Tr 


{{{3a + l)R-3Pk)f' -f) ^dtR, 


?( 2 ) 


(((3/3 + 1)R - 3Pfc) V" - ((3/3 + 2)R - Wk) f + 2/) ^ 


4 °^ 


(3.5) 


Here t = Ink/ko is the dimensionless “renormalization group time”, the primes denote 
derivatives of fkiR) with respect to R and = □ -|- Rfc(n) denotes the regularized kinetic 
term constructed from the spin-dependent coarse-graining operator = —D'^ + E(^gy In 
the nomenclature of [24], the inclusion of non-trivial endomorphisms in the regulator 
function gives rise to a type II regulator. On a spherical background is proportional 
to R and we set 

E(2) = aR , Eyy = j3R. (3.6) 

The traces T(2) and T(o) are with respect to the transverse-traceless and scalar fluctuations, 
respectively. The particular form of the flow equation (3.4) results from either a particular 
choice of the path integral measure in combination with a special choice of regulator for 
the gauge degrees of freedom or, alternatively, through a specific choice of variables on the 
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Spin s 

Eigenvalue A/(d, s) 

Multiplicity Mi{d,s) 


0 

d{d-i) ^ 

{2l+d-l){l+d-2y. 

i\{d-iy. 

1 = 0,1,... 

2 

l{l+d-l)-2 D 

{d-2){d+l)(l-l){d+l){d+2l-l){d+l-3y. 

1 = 2,3,... 

d{d-l) ^ 

2{d-iy.{i+i)\ 


Table 1. Eigenvalues and multiplicities of the Laplacian operator A = —D^ acting on fields with 
spin s on the d-sphere [63, 64]. 


configuration space of the gravitational path integral [60]. For vanishing endomorphism 
the traces coincide with earlier constructions [22, 23]. 

The next step consists in evaluating the traces (3.5). In the case where the flow equation 
is expanded in powers of the background curvature R one typically resorts to the local (or 
early-time) expansion of the heat-kernel which is truncated at the desired order in R. The 
resulting expansion is a good approximation in the regime S> R, but does not give rise to 
the correct IR behavior, <C R, since, by definition, it does not include the non-analytic 
terms reflecting the compact nature of the background. In order to implement construction 
principle 3 to obtain a flow equation with the correct asymptotic for both 3> R and 
< R, we evaluate the functional traces as the sum over the eigenvalues of the differential 
operators. The spectrum A;(d, s) and multiplicities Mi{d, s) of the Laplacian on a spherical 
background are well known and summarized in Tab. I. Introducing the functions IF(j) the 
functional traces can be written as 

OO 

r( 2 ) = Tr [IF( 2 )(A( 2 ) + aR)] = Mi{4, 2) 1^(2) (A,(4, 2) + aR) , 

1=2 
OO 

r(o) = Tr [lT(o) (A(o) + pR)] Mi (4,0) IF(o) (A, (4,0) + /3R) , 

1=0 

where A^^) = —D^ denotes the standard Laplacian acting on fields with spin s. 

In order to evaluate the infinite sums explicitly, it is convenient to switch 
sionless quantities 

r = Rk-\ fk{R)^k^ipk{Rk-^) . 

Moreover, we specify the regulator Rk to be of Litim type [61] 


to dimen- 

(3.8) 

(3.9) 


where 9 is the Heaviside step function. Substituting the traces (3.7) into (3.4) and express¬ 
ing the result in terms of dimensionless quantities then yields 


3847r2 


[ip + Aif- 2r ip') = ^ Mi{4, 2) 


1=2 


di^i ip' 0.2,1 - 2 rip") 

3 ip — ((3a -|- l)r — 3) (p' 


A2) 


+ ^ MK4,0) 


z=o 


gp; ip' a 2 ,i ip" + 03 ,/ + a4,i - 2rip'") 

2ip + {3 - (3/3 -I- 2)r) ip' + ((3/3 -|- l)r - 3)^ (p" 


AO) 


(3.10) 
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Here the dots and primes denote derivatives with respect to t and r, respectively and we 
omitted the arguments of to ease readability. The coefficients in the transverse- 

traceless sum are 

(36 - r{l{l + 3) - 2 + 12a)) , a 2 ,i = ^ (12 - r{l{l + 3) - 2 + 12a - 2)) , (3.11) 

O O 

while the coefficients in the scalar trace are given by 
aM = ^(36-r(/(/ + 3) + 12/3)) , 

a2,i = ^ (432 - 96r(3/3 + 2) + r2 {l\l + Sf - lU/3^)) , 

3"^ (3.12) 

a3,i = -(12-r(/(/ + 3) + 12/3)) , 

04 ,; = ^ (144 - 96r (3/3 + 1) + {f{l + 3)^ - 8l{l + 3) - 48/3(3/3 + 2))) . 

The explicit form of the step functions is given by 

0(2) =0(l-^r (/(/ + 3)-2 + 12a)) , 0(o) = 0 (l - ^ r (/(/+ 3) + 12/3)) . (3.13) 

Eq. (3.10) provides a very explicit illustration how the FRGE (1.1) works: the quantum 
fluctuations are organized with respect to the eigenfunctions of the background Laplacian. 
Each eigenmode thereby gives a specific contribution to the trace whose weight is deter- 
mined through the Hessian T). h From the step functions (3.13) one explicitly sees that 
lowering k (which, for hxed background curvature R, corresponds to increasing r) the eigen- 
modes are integrated out consecutively starting from fluctuations with large values 1: once 
k'^ crosses the corresponding eigenvalues the step-functions ensure that the contribution of 
the associated spherical harmonic drops out from the rhs of the flow equation and does no 
longer contribute to the running of ^pk- 

The structure of the step functions (3.13) also clarifies the role of the endomorphisms 
E(s). The parameters a and (5 allow to adjust the value which the last fluctuation 

mode of the corresponding trace is integrated out. Substituting the /-value of the lowest 
eigenmode and equating the argument of 0/^) to zero yields 

3 1 

^(2),term 2 -\- 3q; ’ ^(0),term ^ • (3-14:) 

For r > r( 5 ) tej-m the corresponding trace vanishes identically. The condition that both 
the spin-two and scalar trace are completely integrated out at the same value rterm can 
be implemented by requiring that the coarse-graining operators 0 ( 5 ) have equal lowest 
eigenvalues (ELE). Equating r( 2 ),term and r(o),term shows that the ELE condition is satisfied 
along the line 

a = /3-^. (3.15) 

In the following, we will implement this condition thereby realizing the guiding principle 
4. This choice is motivated on physical grounds, since it implies that all fluctuations are 
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integrated out at one fixed value of k. It turns out that the ELE-condition can consistently 
be imposed with all other guiding principles. 

In principle the endomorphisms can also be used to integrate out all fiuctuation modes 
on a compact interval by selecting r(2)^term = ^(o),term < oo. For r > 'r(o),term I'hs of the 
flow, encoding the quantum effects, becomes trivial and the flow equation reduces to the 
classical scaling relation. This was the strategy implemented in [54], which constructed 
fixed functions on the compact interval 0 < r < r(o),term = 6- While this compactification 
simplifies the numerical analysis, it will not be implemented here and we will consider fixed 
functions which are well-defined on the entire positive half-axis 0 < r < oo. 


3.2 Smoothing the staircase 

A feature that makes (3.10) particularly difficult to work with is that the PDE possesses 
an infinite number of discontinuities. At the points 

12 , 12 


r2,l = 


I > 2 , 


^0,/ = 


-, /> 0 , 


(3.16) 


12a-h/(l-h3) - 2’ ” ^ 12l3 + l{l + 3)' 

where crosses an eigenvalue of 0(5), the step functions (3.13) remove the corresponding 
eigenmode and the rhs of (3.10) changes discontinuously. This feature is particular to the 
use of the Litim regulator and would be absent in the case of a smooth (e.g. exponential) 
regulator. In order to arrive at a smooth PDE we follow the strategy of [55] and approxi¬ 
mate this staircase-behavior by a smooth function. This is achieved as follows. Firstly, we 
observe that, for r > 0, the sum over eigenvalues is truncated to a finite sum. In order for 
a term to contribute to this series, the corresponding argument of the theta-function has 
to be positive. This gives an r-dependent upper boundary W on the number of terms in 
the finite sum 


12 - r{12a + l{l + 3) - 2) > 0 
for the tensor modes and 

12 -/(3/)r - 12r/3 > 0 ^ 


, ^ V'r((17-48a)r+ 48) - 3r 
I < -^- =: N^2)[r) 


2r 


(3.17) 


I < 


\/3 V-16/3r2 + 3r2 + 16r - 3r 


2r 


=:iV(o)(r) (3.18) 


for the scalar modes. Again, it can be seen that only for r = 0 infinitely many terms are 
contributing. Regarding the 1-dependence, the spectral sum contains only a polynomial in 
I, as can be seen in eq. (3.12). For r > 0 each monomial can be summed using the relation 


M 

E 

71=1 


= 


k 

E 

j=0 


^k—j 

j + 1 


M^'+^, 


2 


(3.19) 


with Bernoulli numbers Bn- Substituting M i—)■ A’(j)(r) yields a smooth approximation 
bounding the staircase-function from above. A smooth approximation that bounds the sum 
from below is obtained by a sum from n = 1 to n = A^^ ~ 1- In practice, we approximate 
the staircase by the average of this two approximations, setting 


Nr 


1 


Nr 


, n I—- 
^ 2 
n=l 


E"‘ + 


\n=l 



(3.20) 
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Applying the summation formula (3.20) to the rhs of the flow equation (3.10) yields 
the desired smooth PDE govering the scale-dependence of 

c\Lp' -I- - 2rip'') 

(3.21) 


(/j + 4(/9 — 2rLp' = 


3(f — ((3a + l)r — 3)(p' 

C\ip' -I- C2<^" -I- C3<y9' -I- C4 ((^" - 2 r^"') 




2ip + {3 — (3/3 + 2)r)ip' + {3j3r + r — 3Yip'' 

The coefficients c* and q are polynomials in r and depend on the endomorphisms 

5(6«r+r’—6)((l8a^+9Q:—2)r^ —18(8o+l)r-|-126) 

~ 691271^ 

~ 5(6Q!r+r—6)((3o+2)r—3)((6 q!— l)r—6) 

*^2 — 691277^ 


(3.22) 


and 


_ ((6/3-l)r-6)(/3(6/3-l)r2+(10-48/3)r+42) 

*^1 ~ 230477^ 

_ ((6/3-l)r-6)(/3(54/32-3/3-l)r3 + (270/32+42/3-35)r2-39(18/3-l-l)r+378) 

*^2 — 460877 ^ 

_ (/3r-l)((6/3-l)r-6)^ 

230477^ 

_ (/3r-l)((6/3-l)r-6)^((9^-l-5)r-9) 

‘'4 460877 ^ ■ 


(3.23) 


The RG flow of the function ^Pkif) is governed by the PDE (3.21). Within this type of 
functional truncations, fixed points are promoted to fixed functions (p*{r) given by globally 
well-defined, fc-stationary solutions of the PDE. Dropping all t-derivatives from (3.21) then 
yields a third-order ODE for the fixed function (p*{r) 

cup' — 2 c 2 rip" Clip' + C 2 ip" — 2cirip"' 


^ip — 2rip' = 


+ 


3(p — {3ar + r — 3)ip' (3/3r -\-r — 3)‘^ip" -|- (3 — (3/3 -(- 2)r)ip' + 2ip 


(3.24) 


with coefficients given in (3.22) and (3.23). Since there is no risk of confusion we will omit 
the asterisk and denote the fixed function only by (p{r). 

One particular effect of the smoothing procedure applied in this subsection is that rhs 
of the flow equation is non-zero on the entire interval r £ [0, oo[ even if the endormorphisms 
are chosen in such a way that the quantum fluctuations in the initial flow equation (3.10) 
are integrated out on a finite r-interval. Thus, implementing the guiding principle 5, we 
insist that a “globally well-defined” solution is required to be regular on the full interval 
r £ [0, oo[. 


3.3 Conformally reduced flow equation 

Based on the full flow equation (3.21), it is straightforward to write down the “conformally 
reduced” approximation where the RG flow is solely driven by fluctuations of the scalar 
mode. Eliminating the contributions from the transverse-traceless fluctuations by setting 
Cl and C 2 to zero yields 


ip + Aip — 2rip' 


Clip' + C 2 ip" + c^ip' -h C 4 {ip" - 2rip'") 

2ipk + (3 - (3/3 -I- 2)r)ip' + (3/3r + r - 3Yip" ’ 


(3.25) 
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with coefficients (3.23). Setting all t-derivatives equal to zero yields the fixed point equation 


4 ^ - 2r^' = _ ci^' + C2ip''-2c,np''' _ 

^ ^ (3/3r + r - 3) V" + (3 - (3/3 + 2)r)ip' + 2ip 


(3.26) 


which is the analogue of the ODE determining the fixed functions in the full system (3.24). 
Comparing the ODE’s of the conformally reduced and full system one finds that they share 
all the characteristic features: both settings give rise to a non-linear ODE of third order. 
Moreover, the coefficient multiplying the highest derivative of is identical, indicating 
that both equations possess the same fixed singularities. On this basis, one may expect 
that both equations also share the same number of globally well-defined solutions. This 
expectation will be confirmed in Appendix A where eq. (3.26) will be analyzed in detail. 


4 Analytical properties of the stationary flow equation 

The ODE (3.21) encoding the fixed functions of QEG in the /(i?)-approximation is a 
complicated non-linear third order differential equation making the construction of exact 
analytic solutions a difficult task. Nevertheless some features of the ODE can be accessed 
by analytic methods. We start with analyzing the fixed singularities in Sect. 4.1 before 
zooming into the asymptotic IR behavior of solutions in Sect. 4.2. 

4.1 Singularity structure 

Being of third order, the ODE (3.21) gives rise to a three-parameter family of solutions 
which exist at least locally. The condition that such a local solution can be extended to 
a global solution, existing for all r > 0, gives rise to constraints on these parameters. If 
the ODE contains fixed singularities, global solutions have to fulfill regularity conditions at 
these points in order to be able to pass through the singularity. This constrains the space 
of solutions and generically one expects that each fixed singularity reduces the number 
of free parameters by one. The existence of isolated, globally defined solutions generally 
requires a total number of three singularities. This condition is met, for example, if the 
ODE has three fixed singularities in its domain of validity. 

In order to determine its fixed singularities, the ODE (3.21) is cast into normal form 
by solving for One then sees that the coefficient C 4 multiplying r(p'"{r) determines 

the fixed singularities of the ODE: zeros of C 4 constitute points where the ODE develops 
a pole. Due to the factor r multiplying one pole is always real and located at 

r = 0. This pole originates from switching to dimensionless quantities and the fact that R 
is not dimensionless. The coefficient C 4 is a forth order polynomial in r, giving rise to four 
additional roots. Thus set of fixed singularities is given by 

9 1 6 

J’sing,! — 0 , 'rsing,2 — ^ ’ ^sing,3 — ^ j ?’sing,4,5 — ^ • (4-1) 

The position of the roots Tsing^, i = 2,3,4, 5 depends on the endomorphism /3 and is 
shown in Fig. 2. Notably, rsing ,2 diverges at /3 = — | and the double pole rsing, 4,5 is shifted 
to infinity if /3 = ^, so that these are distinguished choices for the endomorphism. The 
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Figure 2. Position of the roots rinsing of the coefficient C 4 , eq. (3.23), as a function of the endo¬ 
morphism (3. At the dashed vertical lines one of the roots is located at infinity, implying that the 
number of fixed singularities on the positive half-axis changes at these values of p. 


analysis in Fig. 2 shows that that for /3 = 0 the ODE has two fixed singularities at finite, 
positive values r. By choosing the endomorphism in the interval 0 < /3 < ^ the number of 
hxed singularities is increased to three. For P > \ the number of fixed singularities further 
increases to four. Thus the inclusion of the endomorphism allows to control the singularity 
structure of the ODE to a certain extend. 

Based on the singularity counting argument presented in Sect. 1 . 2 , we then select /? 
such that the number of fixed singularities matches the order of the ODE. In anticipation 
of the numerical analysis, we will impose the condition of ELE and choose 

/3=g> “ = 

In this case, the coefficient C 4 reduces to a third order polynomial which considerably 
simplifies the numerical analysis of the ODE computations. 


4.2 Asymptotic IR behavior of the solutions 

The second property which can be accessed analytically is the scaling of the solutions ^p{r) 
in the IR limit r ^ 00 . In order to zoom into this asymptotic region, we rescale r and (p 
by a constant scale parameter e, 

r = e~°'r, ip{r) = e~^(p{r) , a>0, (4.3) 

and require that the new variables r and p are strictly of the order one. Since a > 0 
the limit e —)• 0 corresponds to r — 00 . Retaining the leading terms in e only, eq. (3.24) 
simplihes to 

1 /. . ^ 2f(p''-(p' 

rp) - f + 6 
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where (p = p>{f) and the primes denote derivatives with respect to the argument. The 


Ihs is independent of the scale parameter a, which can be traced back to the fact that 
the Ihs of eq. (3.24) is invariant under rescalings of the form r i—)• Ar, ip ip. Similarly, 
the rhs is independent of the scale parameter b reflecting that the rhs of eq. (3.24), is 
invariant with respect to the rescaling ip i— \ip, r i—)• r. Notably, the rhs of (4.4) receives 
contributions from the transverse-traceless sector only. The choice /3 = 1/6 reduces the 
order of polynomials Cj so that the contribution of the scalar sector is sub-leading in the 
large-r limit. 

Depending on the ratio of a and 6, eq. (4.4) exhibits three different scaling behaviors; 
classical scaling for b > 2a, quantum scaling for b < 2a and balanced scaling b = 2a. These 
cases will be discussed below. 

a) classical scaling b > 2a 

For b > 2a, the scaling behavior implied by eq. (4.4) is dominated by its classical Ihs 


Aip — 2f(p' = 0 . 


(4.5) 


This equation is solved by 


ipif) = Af"^ 


(4.6) 


with A being a free integration constant. 

b) quantum scaling b <2a 

In this case the asymptotic behavior entailed in (4.4) is dominated by the rhs encoding the 
quantum fluctuations. Disregarding potential singularities caused by the denominator, the 
scaling is governed by the linear equation 


2fip" — ip' = 0 


(4.7) 


which entails that 


3 

(^(f) = Alfa -I- A2 , 


(4.8) 


3 

with Ai and A 2 free integration constants. Notably, the existence of the <p(f) = Aif^ solu¬ 
tion is not intrinsic to the special choice of endomorphisms used to derive (4.4) but appears 
for generic values a and /3 where the scalar sector also contributes to the asymptotic IR 
scaling. 


c) balanced scaling b = 2a 

In this case, the classical and quantum terms balance and eq. (4.4) entails that 



(4.9) 


Solving this non-linear equation exactly is beyond the scope of this work. Instead we check 
if it admits a scaling solution of the form ip{f) = Ar^. Inserting this scaling ansatz into eq. 
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(4.9) yields 


2(2 - a)Af“ 


5a(2o — 3) _2 
8647r2(a + 6)^ 


(4.10) 


For this equation to be satisfied, f has to have the same power on both sides implying 
classical scaling a = 2. For this choice the Ihs vanishes, however, so that o = 2 does not 
constitute to a consistent solution. 

The structure of (4.9) suggests modifying the scaling ansatz including a logarithmic 
term 


(p{f) = Af'^ (1 — r/log(f)) . 


(4.11) 


Inserting this modified ansatz into eq. (4.9) leads to a consistent scaling equation. Expand¬ 
ing the result to linear order in r] then fixes the coefficient in terms of A 


40 

^ “ 552967r2A -h 95 ' 


(4.12) 


Also, in this case the appearance of a logarithmic contribution in the scaling solution 
constitutes a generic feature and is not related to the particular choice of endomorphisms. 

The scaling solutions obtained in the cases a), b), and c) constitute the three admissible 
asymptotic behaviors of potential solutions of the full ODE in the IR. A priori it is unclear 
which phase is realized by a potential fixed function. Nevertheless, this scaling analysis 
provides important information on the IR dynamics of potential fixed functions in QEG. 


5 Numerical construction of the fixed function 


In this section, we use a numerical shooting method to construct solutions of the ODE 
(3.24) which are regular at the fixed singularities. The central result of this construction 
is the unique, isolated and globally well-dehned fixed function shown in Eig. 4. 


5.1 Implementing regnlarity at the fixed singularities 

Our starting point is the ODE (3.24) for the specihc choice of endomorphisms (4.2). The 
analysis (4.1) shows that in this case, the ODE possesses three fixed singularities located 
at finite values r: 

18 , , 
^l,sing — 0 , ^2,sing — ’ ^3,sing — 6 . (^-l) 

Regularity of the solution at these points imposes non-trivial boundary conditions on <p{ri) 
and its derivatives. Casting the ODE into normal form by solving for , and expanding 
the resulting rhs at one of the fixed singularities yields a Laurent expansion of the form 


p"\r) = 


ei{p"{ri),Lp'(ri),ip{ri),ri) 


r - n 


+ regular terms, 


(5.2) 
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with the residues e* given by 


ei 


ip' + Lp 


ip^pp' + \PP'' + ‘^{P')'^ + 


f^V') 


{9p" + ^p' + 2p) , 


62 


1 

16v3'+13</p 


/ 46328 
2535 


PP 


400489 

3510 


PP' + 


14837 

2535 


{P'? - 


+ ^ {Up - 9(/j') {np" - 39p' + imp) , 


2468992 / // , 456000 

32955 P P ^ 28561 



1 / 73 / 

63 = 


f W" 






) “ “ ^p') (l V' - 6 ^?' + p) ■ 

(5.3) 


Here it is implicitly understood that p{r) and its derivatives appearing in e, are evaluated 
at ringing. Keeping p"'(r)\r=ri finite then requires e, = 0. A globally well-defined solution 
has to satisfy the resulting boundary condition (BC) at all three singular points (5.1). On 
this basis, we expect that the BCs fix all three free parameters characterizing the solutions 
of the ODE locally, so that the set of globally well-defined solutions is discrete. 

The next step consists in characterizing the solutions of the ODE by free parameters. 
Imposing that p{r) satisfies the regularity conditions (5.3) ensures that the corresponding 
solution can be expanded in a Taylor series at r^^sing since all its derivatives remain regular. 
In particular a fixed function has an analytic expansion in the UV (r = 0), so that the 
solution locally takes the form 

Pi {r) = go + gir + 9n r” . (5.4) 

n=2 

Substituting this ansatz into the ODE and imposing regularity at ringing by requiring ei = 0 
uniquely fixes the constants gn, n > 2 in terms of {go,gi)- In particular, 

_ 7ffi (I3ffo -F 18gi) -F 7687r^ffo (2gg - 5goffi - 3gf) 

“ 13824772^0(50 + 5i) - 63(35o + 135i) ' ^ ^ 

Thus, implementing the regularity condition ei, we obtain a two-parameter family of (local) 
solutions specified by ( 50 ) 5 i)- 

In a similar fashion, one can perform an analytic expansion of p{r) at the other singular 
points 


P2{r) = ho + bi{r - r2) -h ^ bn{r - r2)"', 

71=2 

(5.6) 

Psir) = Co + ci(r - 5 ^ 3 ) + ^ Cn{r - rg)” . 

(5.7) 


n=2 


Combining the ODE with the corresponding regularity condition again fixes the the coef¬ 
ficients bn and c„ with n > 2 in terms of {ho, hi) and (co,ci), respectively. Eor the c„ this 
solution is again unique, while the occurrence of {p")"^ in 62 makes the solution 62(^01 &i) 
double-valued. In the following we will use the parameterization {go^gi) to characterize 
potential global solutions, which then already satisfy regularity at ringing = 0 . 
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5.2 Extending the local solutions away from r = 0 

The next step consists in extending the two-parameter family of solutions (5.4) to the 
interval r S [0, oo[ implementing the boundary conditions (5.3). In contrast to the example 
discussed in Sect. 2, we do not have an analytic solution of the ODE at our disposal. 
Moreover, the polynomial expansion is difficult to implement to sufficiently high order 
due to the complexity of the ODE. Therefore, we rely on the numerical shooting method 
illustrated at the end of Sect. 2 . Since we need to implement the boundary conditions at 
?’ 2 ,sing and rs^sing, we divide the positive half-axis into three intervals 

h^] 0 ,r 2 [, / 2 =]r 2 ,r 3 [ and/s =]r 3 , oo[. (5.8) 

The shooting algorithm is applied twice: firstly, the solution is constructed on the interval 
Ii- Solutions meeting the boundary condition 62 = 0 are subsequently extended into the 
interval I 2 where the condition 63 = 0 is tested. Finally, solutions satisfying all three 
boundary conditions are numerically extended onto I 3 . Owing to the non-linearity of the 
ODE, a solution might also terminate in a movable singularity so that it does not extend 
up to the hxed singularities marking the upper boundaries of the integration intervals. 

Practically, we implement this search strategy for fixed functions as follows. Because 
of the fixed singularity at rinsing = 0 , the initial conditions for the numerical integration 
need to be imposed at r = ei rather than at r = 0. The initial values for ip{r)\r=ei and its 
derivatives in terms of the free parameters {go,gi) are obtained from the expansion (5.4), 
yielding 

ip{ei) = ipi{ei-,go,gi), ip'{ei) = ip'^{ei-,go,gi), ip"(ei) = go, gi). (5.9) 

We then discretize the (^o, 5 'i)-plane by a lattice and use the initial conditions obtained at 
each lattice point to numerically integrate the ODE from ei to r 2 ,sing — £ 2 - This defines a 
two-parameter family of solutions 

^{r;go,gi), r e [ei, r 2 ,sing - £ 2 ] • (5.10) 

Evaluating the BC 62 , eq. (5.3), based on these solutions provides a map 

{go,gi)e 2 {r 2 - e 2 ;go,gi). (5.11) 

The values for 62 obtained in this way are shown in the top diagram of Fig. 3. Here 
lattice points with 62 > 0 are marked green while points with 62 < 0 are indicated in red. 
The regular points satisfying 62 ~ 0 are given by the black line. This set of regular points 
constitutes a one dimensional subspace that will be referred to as regular line. Every point 
on the regular line can be identified by go, which is the only remaining free parameter. 
Thus implementing regularity of the solution at the hxed singularities rinsing and r 2 ,sing 
reduces the dimension of space of local solutions by two. 

In the next step, we cross the hxed singularity by utilizing that the solutions regular 
at r = r 2 ,sing admit the polynomial expansion (5.6). The two coefficients bo and bi are 
determined by matching the expansion with the numerical solution (5.10) 

y^{r 2 - 62 ] go) = P 2 {r 2 - e 2 ]bo,bi), ip'{r 2 - € 2 ', go) = ¥^ 2(^2 - e 2 -,bo,bi). (5.12) 
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Figure 3. Summary of the results obtained from the numerical shooting algorithm. Solutions 
satisfying the boundary condition 62, characterized by the two free parameters (gojffi) are shown 
in the top diagram. Points with 62(30,51) > 0 are marked green while points with 62(30,31) < 0 
are marked red. The black line indicates the parameters which lead to solutions which are regular 
at rinsing and r2,sing- Along this line 31 is determined as a function of go- The extension of these 
solutions to the interval I2 is depicted in the bottom diagrams. Notably, most of the solution do 
not extend to ranging but terminate in a movable singularity. The point of termination r^ax of the 
numerical solution as a function of go is shown in the left panel. For the solutions extending to 
ra.sing, the regularity condition 63(30(31)) is evaluated in the right panel. The minimum shown by 
63(30) indicates that there is one isolated solution which satisfies all three BCs, 6 ^ = 0 . 


Here we explicitly indicated that the numerical solution depends on the single remaining 
free parameter 30 , stressing that we consider the regular line only. Given the map 

go^{bo,bi) (5.13) 

we evaluate the polynomial expansion at r 2 ,sing + £2 in order to obtain the initial conditions 
for the numerical extension of the regular solutions to the interval I 2 

^{r 2 + £ 2 ) = ‘^ 2 {r 2 + 62 ;3o), 

<^'(^2 + 62 ) = <y? 2(^2 + e2;5o), (5.14) 

g^''{r 2 + 62 ) = (P 2 {r 2 + £ 2 ; go) ■ 

The results obtained from the second numerical integration are shown in the bottom 
diagrams of Fig. 3. Notably, only a small window of points on the regular line, 0.029 < 
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90 

91 

n 

^0 

bi 

Co 

Cl 

full flow 

0.0299 

-0.0364 

0.224 

-0.0154 

-0.0342 

1.3266 

0.7521 

conformally reduced 

0.0528 

-0.0587 

0.152 

-0.0198 

-0.0448 

13.558 

7.1777 

polynomial expansion 

0.0111 

-0.0237 

0.197 

- 

- 

— 

- 


Table 2. Numerical values of the first two coefficients appearing in the polynomial expansions 
( 5 . 4 ), ( 5 . 6 ), and ( 5 . 7 ) of the fixed function at the fixed singularities ( 5 . 1 ). The first line shows 
the coefficients obtained from the full flow equation. The values from the conformally reduced ODE 
( 3 . 26 ) are derived in Appendix A and shown in the second line. For completeness, we also give 
the expansion coefficients of the polynomial /(i?)-expansion obtained in [ 60 ] in the third line. The 
universal product r* = gQ/{ 32 '!rgi) is qualitatively similar in all three cases. 

9 o ^ 0.0315, the solutions can be extended up to raising — £ 3 - All other solutions terminate 
in a movable singularity. Their corresponding end points are displayed in the bottom-left 
diagram of Fig. 3. For the solutions reaching raising — 63 the numerical integration provides 
a map 

go esirs - €3; go) . (5.15) 

The values of e 3 (r 3 — € 3 ; go) are shown in the bottom-right diagram of Fig. 3. The map 
has a pronounced minimum at go ~ 0.0299. At this minimum the BC 63 ~ 0 is satisfied 
within the numerical accuracy of the shooting algorithm. This indicates that there is a 
single, isolated solution which is also regular at r 3 _s ing - 

Finally, the isolated solution is extended to the interval I3. Following the strategy 
employed when crossing the singularity at r 2 ,sing) the numerical solution is matched to the 
polynomial expansion (5.7) 

9^ir3 - 63; go) = ip3ir3 - e3;co,ci), (^'(rs - 63150 ) = ¥^ 3(^3 - e3;co,ci). (5.16) 

This determines the values of cq and ci for the regular solution. Following the strategy 
implemented in eq. (5.14), the parameters co,ci are converted into initial conditions for 
a numerical integration on I3. The solution <f*{r) is then completed by carrying out this 
numerical integration and matching to the scaling regimes discussed in Sect. 4.2.^ The 
numerical values of all expansion coefficients for the regular solution are summarized in 
the first line of Tab. 2. The resulting isolated and globally well-defined fixed function (/5*(r) 
is shown in Fig. 4 and its properties are discussed further in Sect. 5.3. This solution is 
the first proper fixed function obtained within four-dimensional Quantum Einstein Gravity 
and constitutes the central result of this work. 

We close this subsection with some technical comments on the implementation of the 
shooting algorithm. Notably, the position and width of the plateau displayed in the bottom- 
left diagram of Fig. 3 has a mild dependence on the parameters entering the numerical 

^Based on the values co, ci, the the solution on the interval I 3 can also be obtained by recursively solving 
the polynomial expansion (5.7). Carrying out the expansion up to A = 10, one finds that this expansion 
agrees with the numerical solution within its radius of convergence. 
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Figure 4. The isolated, globally well-defined fixed function (/?*(?') constructed from imposing the 
three regularity conditions ( 5 . 3 ). The fixed function is bounded from below, comes with a positive 
cosmological constant A* and Newton’s constant 5*, and displays a global minimum at rmin ~ 2 . 5 . 


integration. It is conceivable that it reduces to a single point (spike) when implementing a 
higher numerical precision. Owing to the fact that we are using a polynomial expansion at 
the singular points to continue the solutions through a fixed singularity the entire construc¬ 
tion is very robust under variations of the auxiliary parameters ei, €2 and 63 . Changes in 
Cj induce small changes (within error-bars) of qq that propagate into the coefficients 6 o)^i 
and Co, Cl. These quantitative changes do not affect the qualitative picture shown in Fig. 
3, however. 

5.3 Properties of the complete fixed function 

The proper hxed function (/?*(r) satisfying all BCs (5.3) is shown in Fig. 4. Here, we will 
summarize its main properties. 

Asymptotic expansion in the UV 

At r = 0, which corresponds to fc —)■ 00 if the background curvature is held fixed, the fixed 
function admits the analytic expansion (5.4). The first two expansion coefficients are given 
in the first line of Tab. 2. In order to facilitate the comparison with results obtained in 
finite-dimensional approximations, it is convenient to compute the corresponding values 
of the dimensionless Newtons constant g* = —{lQTTgi)~^ and cosmological constant A* = 
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—g{)/{2gi) together with the universal (gauge-independent) product r* = g^ A* 

A* = 0.411, 5* = 0.547, n = 0.224. (5.17) 

These values are in good agreement with the results obtained from expansions of flow 
equations truncated at finite order [8, 16, 17, 65]. In particular r* fits well to results 
obtained within “geometric” flow equations [31, 60]. We take this agreement as strong ev¬ 
idence that, despite the different boundary conditions imposed on the flow equations, the 
polynomial expansions and the fixed function (p*{r) actually describe the same universality 
class. 

Consequences of the global minimum 

The fixed function possesses a global minimum located at Cmin ~ 2.500. At this minimum 
<7’*(?'min) < 0. This feature together with the resulting non-monotonicity distinguishes our 
solution from earlier constructions [55, 56], where turned out to be positive definite 

and monotonically increasing. 

As a consequence our solution actually passes the redundancy test [58], requiring that 

E 4 {r) = 2(y9*(r) - r ^p'^{r) (5.18) 

has a zero for some value r. From Fig. 4 we then deduce 

.^ 4 ( 0 ) = 2^0 > 0, F'4(?’min) = 2(/7*(rmin) < 0. (5.19) 

Since E 4 ^{r) is continuous this implies that Eii[r) has a zero within the interval r S [0, Tmin]- 
This provides a strong indication that linear perturbations of the fixed function cannot be 
absorbed by a field redefinition and constitute genuine eigenperturbations of the fixed func¬ 
tion. 

Asymptotic behavior in the IR 

Finally, it is interesting to investigate the asymptotic behavior of </?*(r) for r —)• oo. In 
order to determine the scaling regime, we first approximate the large-r behavior through 
the ansatz 


ip{r) = Ar"^ {l + uir ^+U 2 r , (5.20) 

neglecting further, sub-leading terms. The leading behavior is motivated by the clas¬ 
sical and balanced scaling asymptotics encountered in Sect. 4.2. Fitting eq. (5.20) to the 
numerical solution obtained on the interval r £ [6,9] yields the parameters 

= 0.07705 ± 0.00032 , 

uf = - 2.07514 ± 0.05399 , (5.21) 

uf = - 6.36855 ±0.25897. 

Fig. 5 shows that the ansatz indeed provides a good approximation of the numerical solution 
in this regime. 
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Figure 5. The numerical solution (red points) and the non-linear fit of the ansatz (5.20). The 
values of the fit parameters are given in eq. (5.21). 


At this stage it is also illuminating to substitute the ansatz (5.20) into the ODE (3.24) 
and determine the coefficients A, ui and U 2 from a large-r expansion. Surprisingly, the 
ansatz provides an asymptotic solution if and only if the endomorphisms satisfy 

43 -h 300a -h 360a^ - 36/3 -h 108/3^ = 0 . (5.22) 

This condition is not satisfied for the choice (4.2) and becomes compatible with the ELE 
condition (3.15) for two specific values /3± = (18 ± \/285)/78 only. Therefore we conclude 
that the asymptotic behavior of ^p*{r) is described by the scaling solution obtained in the 
balanced regime and involves also non-analytic terms not present in the ansatz (5.20). This 
is actually in agreement with earlier works [56, 58, 66] which also observed the presence of 
non-analytic terms in this asymptotic scaling regime. 

6 Conclusions and discussion 

The existence of a NGFP on the gravitational theory space is a crucial ingredient for con¬ 
structing a consistent and predictive quantum theory of gravity along the lines of Wein¬ 
berg’s Asymptotic Safety scenario. In this work, we used the functional renormalization 
group equation for the effective average action Tfc, (1.1), to establish the existence of a 
suitable fixed function in the realm of /(i?)-gravity in four spacetime dimensions. De¬ 
spite the vast body of earlier works [55, 56, 58, 67], this is the first time that a suitable 
fixed function has been found and its construction constitutes an important step towards 
generalizing finite-dimensional truncations of the gravitational RG flow to approximations 
containing an infinite number of scale-dependent coupling constants. Our solution is glob¬ 
ally well-defined on the positive real axis and its UV-expansion at r = 0 yields a positive 
cosmological constant and Newton?s constant. The structure and properties of the fixed 
function is in good agreement with the fixed points solutions found at the level of poly¬ 
nomial /(ii)-truncations [60]. This construction constitutes an important step towards 
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generalizing finite-dimensional truncations of the gravitational RG flow to approximations 
containing an infinite number of scale-dependent coupling constants. 

The fixed function is based on non-linear ODE of third order, which construction prin¬ 
ciples discussed in the introduction. Owing to the inclusion of non-trivial endomorphisms 
in the regulator the ODE possesses three fixed singularities at finite r = R/k"^. This fea¬ 
ture restricts the space of globally well-defined solutions to one isolated fixed function. The 
condition of “equal lowest eigenvalues” thereby plays an important role for the existence 
of this solution: numerical analysis indicates that once the transverse-traceless and scalar 
sectors of the flow equation are out of tune, the numerical solutions typically terminate in a 
movable singularity. Combining a scaling analysis of the ODE together with the numerical 
solution shows that /*(i?) = AR?‘{1 + r^logii -|- ...), matching the IR behavior of earlier 
constructions [55, 56, 58, 67]. 

In Appendix A we analyzed conformally reduced version of the full flow equation. 
The resulting ODE shares many properties with the full flow. In particular it is of the 
same order and exhibits the same fixed singularities as the full system. At the same 
time the conformally reduced system is significantly simpler which allows investigating the 
existence of fixed functions by employing complementary techniques. As the central result 
we established that the conformally reduced equation also exhibits one globally well-defined 
isolated fixed function whose properties are qualitatively similar to the ones found in the full 
system (cf. Eigs. 4 and 7). The result establishes that conformally reduced approximation 
is a very useful setting when searching for fixed functions of the gravitational flow equation. 

At this stage, we have not yet analyzed the structure of admissible, globally well- 
defined, inhnitesimal deformations of the fixed function. Substituting the corresponding 
ansatz describing these deformations into the PDE (3.21) leads to a linear third-order 
ODE which also contains the (possibly complex) scaling exponent of the deformation. 
In principle the corresponding free parameters can again be determined by applying the 
combination of analytical and numerical techniques employed in the construction of the 
fixed function. This numerically very demanding, challenge is, however, beyond the scope 
of the present work. The match between the fixed function and its finite-dimensional 
sibling let us expect that the full functional also shares the stability properties of the 
polynomial approximation. Presupposing that all relevant eigenvalues can be extended to 
global deformations, this implies that there should be three relevant deformations. This 
expectation that there is a finite number of relevant directions is also supported by the 
general arguments [57]. As an important first check in this direction, we explicitly verified 
that our solution passes the redundancy test [58] , indicating that the fixed functional might 
indeed have non-trivial deformations which cannot be absorbed by a field redefinition. 

Let us again stress that carrying out the extension of the fixed points seen in finite¬ 
dimensional truncations to the level of fixed functions is an important step in the Asymp¬ 
totic Safety program: the functional approximations are much more sensitive to structural 
details of the theory, as, e.g., the functional measure used to construct the flow equations. 
These details have a crucial influence on the singularity structure of the ODE describing 
the fixed function and may provide important insights on conceptual properties required to 
achieve Asymptotic Safety. In this context, we find it encouraging that the geometrically 
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motivated flow equation actually admits an interesting fixed function solution. 

Notably, our findings are not in tension with the recent results [59] where no fixed 
function for four-dimensionally conformally reduced gravity has been found. The flow 
equation and search strategy implemented in this work differs from our setting in var¬ 
ious ways, which can be made explicit based on the construction principles detailed in 
the introduction. Firstly, the flow equation constructed in [59] utilizes a conformally flat 
background while our equation has been derived for a compact background with positive 
curvature. Secondly, the analysis in [59] implemented the modified split ward identities of 
the effective average action to reduce the hxed point equation to a linear PDF. While this 
certainly simplifies the corresponding fixed function analysis the special choice of regulator 
or the implementation of split symmetry restoration at the level of the truncation may be 
overconstraining the system such that no globally well-defined solution may be found. 

In this light, it would be very interesting to investigate to which extend the construction 
principles of our proposal are strictly necessary for the existence of a satisfactory fixed 
function. In particular it would be interesting to investigate if a flow equation based on 
the exponential parametrization of the metric fluctuations 

Q^iu — dfip ^ ’ (6-^) 

recently advocated in [68-70] also gives rise to hxed functions at the level of the f{R)- 
approximation. Compared to previous constructions, this parameterization corresponds 
to a change in the path integral measure. Combined with a clever choice of gauge-hxing 
this may lead to a second order ODE describing the hxed function. This structure is much 
simpler to analyze since hxed functions may be identihed based on a single shooting instead 
of the double-shooting method required in this work. In [69, 70] it was already observed 
that this choice signihcantly improves the IR behavior of the gravitational RG how. 

Acknowledgements We thank M. Reuter, A. Bonanno, T.R. Morris and K. Falls for inter¬ 
esting discussions. The work by F.S. is supported by FOM grants 13PR3137 and 13VP12. 

A Fixed functions of the conformally reduced flow equation 

In this appendix, we investigate the hxed functions entailed by the conformally reduced 
how equation (3.26). We establish that the conformally reduced setting gives rise to one 
hxed function whose properties are qualitatively very similar to the one constructed from 
the full how equation. The main advantage of the conformally reduced system is that it 
is structurally signihcantly simpler, so that the structure and properties of the resulting 
solution can be corroborated using both analytical and numerical methods. Following 
our choice of endomorphisms from the main part of the paper, we will work with jd = \ 
throughout this section. 

A.l Fixed functions from the shooting method 

We start by analyzing the singularity structure of (3.26). Since the positions of the hxed 
singularities are determined by the scalar trace, the discussion of Sect. 4 can be applied to 
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(3.26) and the the fixed singularities are again given by (5.1). The boundary conditions 
that fixed functions have to satisfy are determined by Laurent expanding according to (5.2) 
and explicitly read 


7o9" lip' 256 2 / // / N 

ei = — + -^-+ 3(/:3 + 2ip) , 


9 


96 


97 


,, .... 2567r2 (13v9 - 9ip') il2ip" - 39ip' + 169vj) 

go —- ip — - p -j- - 

65^ 90^ 7605 

eg = 9 {if' - \2ip") - 5127r2 {ip - 3^9') (iS^j" - ^ip' + ip) . 


(A.l) 


Any regular solution satisfying the BCs (A.l) can be expanded in a Taylor series at 
^’i.sing according to (5.4), (5.6) and (5.7). Substituting these expansions into (3.26) allows 
for determining all coefficients bn and Cn with n > 2 in terms of the Hrst two coefficients 
{go,gi), ibo,bi) and (co,ci), respectively. 

In order to construct globally well-dehned solutions of (3.26) we repeat the numerical 
shooting algorithm used in Sect. 5.2. Beginning at r = 0 we have a two-parameter family 
of solutions characterized by {go, gi). To extend these solutions to [0, oo[, we again split the 
positive half-axis into the intervals (5.8). The first shooting on Ii is realized by numerically 
integrating from r = ei to r’ 2 ,sing — ^2 using the initial values (5.9). After discretizing the 
(( 70 j 5 i)-plane we perform for every lattice point the integration and evaluate the BC 62 
with it. This provides a map 


{90,91)e2{r2 - e2; go, 9 i) ■ (A.2) 

The resulting values 62 are visualized in Fig. 8 , where lattice points with 62 > 0 are 
marked green while points with 62 < 0 are tagged red. Regular points have to satisfy 
62 ~ 0. Clearly, there is a line of points separating two regions of different signs (i.e. 
different colors in the diagram) giving rise to a line of regular solutions. The points on this 
line are again parametrized by go- 

Next, we repeat the shooting procedure on I 2 . Utilizing the polynomial expansion 
at r’ 2 ,sing) initial values for the numerical integration are again determined by matching 
according to (5.12). Integrating from r 2 ,sing + £2 to r-g^ging — eg again reveals a window of 
5 o-values where the corresponding solutions extend up to rg^ging. This regular window is 
qualitatively similar to the one displayed in Fig. 3 and located at 0.052 < 50 ^ 0.054. The 
unique value for go minimizing the last BC eg is listed in Tab. 2 and the corresponding 
globally well-defined fixed function ip*{r) is displayed in Fig. 7. It shares all features of 
the result obtained from the full flow (3.24), including a minimum located at rmm ~ 2.153. 
This result underlines the importance of the conformal mode for the existence of fixed 
functions. 

A.2 Fixed functions from bootstrapping method 

We close this section by presenting an alternative bootstrapping strategy relying on poly¬ 
nomial expansions only [71]. In the conformally reduced setup all expansion coefficients gn, 
bn and Cn are much smaller in size enabling us to go to much higher orders in the Taylor 
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Figure 6 . Solutions regular at r = 0 are parametrized by go and gi. Evaluating these solutions at 
the BC 62 , eq. (A.l), yields the data points depicted above. Points with 62 > 0 are marked green 
while points with 62 < 0 are tagged red. There is a line of points separating two regions of different 
signs (i.e. different colors) giving rise to a line of regular solutions satisfying 62 ~ 0. The black line 
shows gi as a function of go obtained by imposing the bootstrap condition ( 7 jv+i(ffo, ffi) = 0 for 
N = 17 and solving for gi. The two approaches show very good agreement. 


expansions, especially at r = 0. Thus, we start by truncating the Taylor series (5.4) to a 
finite polynomial of the order N 

N 

7^i{r;go,gi) = go + gir + ( 50 , 51 )?’’^ ■ (A.3) 

n=2 

If the Taylor series (5.4) has a finite radius of convergence Ri and if additionally r 2 ,ging < 
i?i, then (A.3) can be substituted into the BC 62 . Again, this provides a map of the 
form (A.2). The condition e 2 {go,gi) = 0 relates go and gi analytically. Unfortunately, 
this relation grows quickly in complexity with increasing order N and thus, is of no real 
practical use. However, if the finite Taylor polynomial (A.3) serves a good approximation 
at the order N, the next coefficient should by very small, i.e. gN+i{go, gi) ~ 0. Thus, a 
much simpler condition that can be imposed is given by 

gN+i{go,gi) = 0, (A.4) 

which constitutes the central idea of the bootstrapping strategy (or “minimal procedure” 
in the terminology of [71]). Being a rational function in go and gi, at every order new 
solutions of (A.4) emerge, which are mostly complex or singular and thus, discarded as 
spurious. 

Pushing the order up to = 17, only one real and non-singular solution of g'lg = 0 has 
been found and is displayed in Fig. 6 (black line). Most remarkably, the analytic solution 
go gi is in perfect agreement with the regular line obtained from the shooting algorithm. 
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Figure 7. The isolated, globally well-defined fixed function of the conformally reduced system 
constructed from imposing the three BCs (A.l). The fixed function is bounded from below, comes 
with a positive cosmological constant A* and Newton’s constant g,, and displays a global minimum 
at Tmin Ri 2.153. 


Moreover, no information from the second pole has been incorporated yet. This can be 
done by repeating the bootstrap strategy at r’ 2 ,ging and imposing 

bM+i{bo,bi) = 0 , (A.5) 

where M is the order to which the expansion (5.6) has been truncated. Imposing (A.5) 
yields a line bo i—;■ 6i in the {bo, 6i)-plane (spurious solutions are discarded), c.f. Fig. 8 (black 
line). In contrast to the expansion at r = 0, no convergence of the line can be observed 
revealing the asymptotic nature of the expansion at r 2 ,sing- Thus, we keep the order M 
low (M = 3 at most) and “transport” the previously obtained regular line from r = 0 to 
^ 2 ,sing into the (6o) ^i)-plane. This is done by substituting the analytic solution go e-j- gi 
into the Taylor polynomial (A.3) and simultaneously setting r = r 2 ,sing- The transported 
line is defined as 

go ^ {g^i{r2,smg-,go,gi{go),g^i{r2,smg-,go,gi{9o)) = {bo{go),bi{go)), (A.6) 

and shown (blue line) in Fig. 8. Inspecting Fig. 8 shows that the transported line (blue 
line) and the solution of (A.5) (black line) intersect at 

bo « -0.0185 6i « -0.0379, (A.7) 
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Figure 8. The black line shows hi as a function of 6o obtained by imposing the bootstrap condition 
(A.5), i.e. bM+iibo, bi) = 0, for M = 3. Due to the asymptotic nature of the polynomial expansion 
(5.6) at r = r 2 ,sing, it is impossible to go to higher orders. The blue line is the result of the bootstrap 
condition gN^i{gQ, gi) = 0 ior N = 17 transported to r 2 ,sing by using the Taylor expansion. The two 
lines intersect at bo Ri —0.0185 and bi fv —0.0379 corresponding to go fv 0.0444 and gi fv —0.0512. 
The bootstrapping strategy confirms the numerical results and thus, the existence of globally well- 
defined fixed functions for the conformally reduced system. 


corresponding to 


go « 0.0444, gi « -0.0512 . 


(A.8) 


Comparing with the values from the shooting method in Tab. 2 shows a good qualitative 
agreement. In summary, the bootstrap strategy provides a further argument in favor of the 
existence of a globally well-defined fixed function. However, due to the asymptotic nature 
of the expansion (5.6) a precise determination of the coefficients characterizing the fixed 
function is impossible. Thus the shooting method employed in the main part of this work 
constitutes the most reliable method for constructing fixed functions in this setting. 
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